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Abstract 

An important challenge in theoretical ecology is to find good-coarse grained repre- 
sentations of complex food webs. Here we use the approach of generalized modeling 
to show that it may be possible to formulate a coarse-graining algorithm that con- 
serves the local dynamics of the model exactly. We show examples of food webs with 
a different number of species that have exactly identical local bifurcation diagrams. 
Based on these observations, we formulate a conjecture governing which popula- 
tions of complex food webs can be grouped together into a single variable without 
changing the local dynamics. As an illustration we use this conjecture to show that 
chaotic regions generically exist in the parameter space of a class of food webs with 
more than three trophic levels. While our conjecture is at present only applicable 
to relatively special cases we believe that its applicability could be greatly extended 
if a more sophisticated mapping of parameters were used in the model reduction. 



1 Introduction 



Over the recent decades there has been significant progress in ecological mod- 
eling. However, this progress manifests itself mostly in the description of small 
systems containing only few species. By contrast, to predict the dynamics of 
large ecosystems remains an important open challenge. In the context of ma- 
rine ecosystems it has often been pointed out that several obstacles have to be 
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overcome in order to model even the planktonic food web, which forms their 
backbone, realistically . 

One obstacle that is encountered in plankton modeling is the sheer diversity 
of marine life. The traditional approach of ecological modeling is to describe 
every single population by at least one differential equation. If we applied 
this approach to a marine food web, we would obtain a huge system, which 
would be prohibitively difficult to study both analytically and numerically. A 
promising alternative is therefore to model the food web not on the level of 
the population, but on a coarse-grained level at which every variable describes 
several similar populations. 

Coarse-grained descriptions have a long tradition in plankton modeling. The 
earliest models, the N-P-Z food chains, consisted only of three equations de- 
scribing, nutrients N, phytoplankton P, and zooplankton Z (see e.g. [H). It 
is clear that models at this level of abstraction cannot capture the dynamics 
of the planktonic food web in any detail. However, they have the advantage 
that the categories N, P, and Z can be clearly defined. Many current models 
use more refined categories such as guilds and functional groups. In this case 
it can be unclear if a given population should belong to one category, or the 
other, or should be a category on its own [2|. The traditional solution is to 
group populations that share certain biological characteristics such as feed- 
ing behavior, metabolism, or activity cycle. Some newer approaches use the 
tools of graph theory to identify cliques or communities that hold a similar 
topological position in the food web [jjj]. The former approach ensures that 
populations that are grouped together in one category, fill a similar role in the 
ecosystem; the latter, that they hold a similar topological position. However, 
does either of the one imply that the grouped species have a similar impact 
on the dynamics? 

In the formulation of a coarse-grained food web model the goal should be 
to group species into categories such that the dynamics of the system is not 
changed by the coarse-graining, at least, not qualitatively. However, to check 
whether this condition is satisfied does not seem to be feasible as it requires 
knowledge of the dynamics of the original system which we set out to determine 
in the first place. Nevertheless, if an approach is available to extract at least 
certain dynamical properties from the full model we can check whether these 
dynamical properties are conserved in the coarse-grained model. 

In this paper paper we apply the approach of generalized modeling to ex- 
plore the effect of coarse-graining on classes of simple food webs. Generalized 
modeling can extract certain dynamical properties, i.e., the local dynamics 
around steady states, very efficiently. One advantage of generalized models is 
that a single generalized model does not describe a single food web, but the 
whole class of food webs which share a similar topology. This enables us to 
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check whether a given coarse-graining step qualitatively preserves the local 
dynamics in a given class of food webs. Based on these results we conjecture 
that certain food webs can be coarse-grained to food chains with the same 
number of trophic levels. This allows us to generalize results on food chains to 
a class of food webs. In particular we present evidence that parameter regions 
in which chaotic dynamics take place, should exist in many food-chains with 
more than three trophic levels. 

The paper is organized as follows: We start by considering some simple exam- 
ples in Sec. 12.11 The insights gained in this section lead to the formulation of 
some conjecture which is stated in Sec. EJ This conjecture is then applied to 
generalize previous results on chaotic parameter regions in food chains to a 
class of food webs. Finally, we summarize and discuss our findings in Sec. HI 



2 Two illustrative examples 



Let us start by using generalized models to study two simple examples. The 
first example is a generalized predator-prey system with one predator. We use 
it to give a brief introduction to the approach of generalized modeling. There- 
after we study a system of two populations of predators competing for one 
prey population. In this second example we discuss the generalized modeling 
only briefly, but focus instead on the results of the modeling process: In certain 
cases the local dynamics of the system is identical to the one observed in the 
predator-prey system. So, in order to capture the dynamics of the system it 
is sufficient to describe the two populations of predators by a single variable. 



2. 1 The generalized predator-prey system 



Let us study the dynamics of a generalized predator prey system. The system 
consists of a prey population X and a predator population Y. In the absence 
of the predator, the prey population grows at the rate S(X). If the predator is 
introduced it consumes prey at the rate G(X, Y) and loses individuals because 
of natural mortality at the rate M{Y). In the following we will not restrict S, 
G and M to specific functional forms, but study the generalized model given 
by 

X = S(X)-G(X,Y) 
Y = G(X, Y) — M{Y) 
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where S, G, and M are arbitrary positive functions. The equations seem to 
suggest that all the biomass that is consumed by the predator is actually con- 
verted to predator biomass. Note however, that any non-vanishing conversion 
efficiency that may exist in a given system can always be set to unity by 
choosing the units of predator or prey biomass appropriately. 

In conventional modeling the first step of model analysis is often to compute 
the steady states of the system under consideration. In the generalized model 
this is impossible with the desired degree of generality. But, there is in fact no 
need to compute steady states: The generalized model describes a whole class 
of systems, which covers the state space densely with steady states. In other 
words, for every steady (X*,Y*) we observe in nature there is at least one 
specific model of the form of Eq. ([I]) that has a steady state exactly at this 
point. In fact there is even a whole family of infinitely many specific models 
that have a steady state at (X*,Y*). 

Let us assume we are interested in an arbitrary steady state (X*,Y*). We 
start our investigations by normalizing the state variables and the biomass 
fluxes to unity. For this purpose we define new variables 

X Y 

x = ™» y = v* ^ 



and functions 

*)=®,^)=^,^)=« (3) 

where S*, G*, and M* denote the fluxes in the steady state S(X*), G(X*, Y*), 
and MiY*). 



Substituting the definitions into Eq. (JTJ yields 

(4) 



x = frs(x) - ^g(x,y) 



y = y^9(x,y) - Y^ m (y) 



The structure of this system of equations is still the same as in Eq. (JTJ), but 
some unsightly pre-factors have appeared. To get rid of these we define the 
constants 



S*_ _ (G*_ _ G* _ M* 



In order to verify the equalities on the right-hand side of these equations 
substitute the normalized steady state (x*,y*) = (1,1) into Eqs. (j4j), use 
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s(l) = g(l, 1) = m(l) = 1, and demand for consistency x = y = 0. 
Using a x and a y we can write Eqs. (0J as 

x = a x (s(x) - g(x,y)) 

(6) 

y = a y (g(x,y) - m(y)) 

In these equations a x and appear as parameters. Indeed they can be treated 
exactly like the parameters that are commonly used in conventional models. 
And, just like them they have an intuitive interpretation: From the definition, 
Eq. (0), we can see that they denote the per capita gain and loss rates in the 
steady state. In other words, a x is the turnover rate of prey, while a y is the 
turnover rate of the predator. 

In order to investigate the local dynamics around the steady state we need to 
compute the system's Jacobian, which constitutes a local linearization of the 
system around the steady state under consideration. In case of our predator- 
prey system the Jacobian can be written as 

tt;j.(S X (^ X ) OL X (jy 

J 

v a y g x a y (g y - m y ) 




where roman indices indicate derivatives of the respective function in the 
steady state. Specifically 



s x = ds(x)/dx\ 1 , (8) 
g x = dg(x,y)/dx\ {11) , 
g y = dg(x,y)/dy\ (hl) , 
m y = dm(y) / dy\ 1 . 

Like the alphas s x , g x , g y and m y are constants for a given steady state and 
can therefore be considered as parameters. In order to understand how these 
parameters are interpreted let us consider how the normalization procedure 
would affect a specific function. Suppose one of the functions in the model, say 
M(Y), was linear. In this case the corresponding normalized function were a 
linear function with slope 1, regardless of the slope of M(Y). Hence, for every 
linear M(Y) the corresponding parameter would be m y = 1. In fact, for an 
arbitrary monomial M(Y) = AY P the corresponding parameter is m y = p. 
For a general function M(Y) we can interpret the corresponding parameter 
as a measure of sensitivity of the fluxes to the variables |4j . 
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In ecology the sensitivity of top-predator mortality to top-predator abundance 
m y is a well known quantity: the so-called exponent of closure. Likewise, the 
other parameters that appear in the predator-prey system have a well-defined 
ecological meaning. The parameter s x denotes the sensitivity of the popula- 
tion growth of the prey with respect to prey abundance. If the prey is not 
limited by other external factors apart from predation the dependence should 
be linear (s x = 1). If by contrast the prey population is strongly limited by ex- 
ternal factors its growth rate should be insensitive to the abundance (s x = 0). 
Perhaps the most important parameter for the dynamics is g x , which denotes 
the dependence of predation on prey abundance. A Lotka-Volterra (i.e., mass 
action) model corresponds to g x = 1. If predator saturation is taken into ac- 
count then the predation rate is less sensitive to prey abundance, and we find 
< g x < 1 depending on the strength of the saturation. Finally, the parameter 
g y describes the dependence of the predation rate on predator abundance. In 
most models this dependence is assumed to be linear (g y = 1). 

Let us now return to the stability of steady states. The Jacobian governs the 
time evolution of the system close to the steady state. If the real parts of all 
eigenvalues of the Jacobian are negative then trajectories that start close to the 
steady state converge to it exponentially. In this case the equilibrium is said to 
be locally asymptotically stable. If parameters are changed the stability may 
be lost if the variation causes eigenvalues to acquire positive real parts. Since 
the Jacobian is a real matrix this can happen in either of two different ways: A 
single real eigenvalue crosses the imaginary axis in the origin of the complex 
plane, or two eigenvalues cross the imaginary axis as a complex conjugate 
pair. The first case corresponds to a bifurcation of saddle-node-type, such 
as fold, pitchfork, and transcritical bifurcations, which in general change the 
number and/or stability of steady states. The second case corresponds to a 
Hopf bifurcation, which, at least transiently, gives rise to oscillatory behavior. 

From the Jacobian, Eq. ([7]), we find that in our model bifurcations of saddle- 
node-type occur if 

g x = s x - (9) 

m y 



A Hopf bifurcation occurs at 



01 

g* = s x -(g y -m y ) (10) 



if (g x — s x )m y + s x g y > 0. Note that only the ratio r = ^ appears in the 
bifurcation conditions. This is reasonable as the absolute values of the turnover 
rates can be changed by means of a timescale renormalization. 
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1.0 1.0 



Fig. 1. Bifurcation diagram of a predator-prey system, depending on the sensitivity 
of the predator to prey abundance g x , the timescale separation r, and the exponent 
of closure m y . Steady states are stable in the top- most volume of parameter space. 
Stability is lost in a Hopf bifurcation (red surface) or in a bifurcation of saddle-n- 
ode-type (blue surface). The sketch in the top-right corner indicates the topology 
of trophic interactions (Parameters: s x = 0.5, g y = 1.) 

The critical parameter values at which the bifurcations occur form hyper- 
surfaces in parameter space. These surfaces are shown in Fig. HJ Every point 
in this diagram corresponds to steady states in a class of models with similar 
local dynamics. The steady states that fall into the topmost volume of param- 
eter space are stable; all others are unstable. Destabilization can occur in a 
Hopf bifurcation (red) or in a bifurcation of saddle-node- type (blue). Among 
other things the diagram shows that a high sensitivity of predation to prey 
abundance (high g x ) has a stabilizing effect on the system. 

Before we move on to the next section, let us recapitulate what we have 
achieved: Starting from a generalized model, Eq. ([1]), we have parameterized 
the Jacobian (and therefore the local dynamics) around every possible steady 
state in the class of models under consideration. Apart from the analytical 
computation of bifurcation points, which we have used here, the Jacobian 
can be analyzed in a number of ways: For instance we can use Monte-Carlo 
sampling to identify the parameters to which stability is most sensitive or 
we can use an optimization approach to find the most stable parameter set 
0; 0] . However, in this paper, we will show that certain topologically different 
food webs can yield exactly the same bifurcation diagrams, which is best done 
analytically. 

2.2 A generalized model of competing predators 

Let us now consider a system in which two populations of predators feed on 
the same prey population. Using essentially the same notation as above, a 
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generalized model for the system can be written as 



X = S(X) - Gi(X, Y x ) - G 2 (X, Y 2 ) 

Yi = Gi(X,Yi) — M X (Y X ) (11) 
Y 2 = G 2 (X,Y 2 )-M 2 (Y 2 ) 



The analytical treatment is very similar to the previous example, for this 
reason we skip some details of the normalization and jump directly to the 
equation corresponding to Eq. (jlj), which now reads 

x = j£s(x) - ^-g x (x,y x ) - ^fg 2 (x,y 2 ) 

y^^g^y^-Mfm^y,) (12) 

V2 = ^92(x,y 2 ) - y£m 2 (y 2 ) 



In analogy to the predator prey system we can define 

J-l J-l 12 J-2 



Note that in the first equation of Eq. ffl2|) there are now three terms. We 
therefore, define 

a r = = 14 

X* X* X* 



and 



G x 



G x * 



Q * ' P 



G 2 * 
G x * + G 2 * 



(15) 



which allows us to write Eq. (fl"2l as 

x = a x (s(X) -f3g x {x,y x ) - (3g 2 (x,y 2 )) 

y\ = oi yl (g x (x,y x ) - m x (y x )) (16) 
V2 = oc y2 (g 2 (x,y 2 ) -m 2 {y 2 )) 
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In analogy to the predator-prey system we compute the Jacobian 



a x (s x - /3g ljX - f3g 2 , x ) -a x [3g hy -a x pg 2 , y ^ 

y a y 2g2,x %2(^2,y ~ ^2,y) j 



(17) 



where now a^f, = Note that a steady state in which all populations have 
positive size does generally exist in the system, as the canonical scenario of 
competitive-exclusion only occurs for g 1>y = m l y and g 2y = m 2,y 0]- 

The condition for bifurcations of saddle-node-type now reads 



)(92 

+P92, x m 2:y (g lty - m lty ) = 



(18) 



Our aim is to compare this condition to the corresponding condition in the 
predator-prey system. This comparison is complicated by the fact that not 
only the state spaces but also the parameter spaces of the two models dif- 
fer in dimensionality. We therefore have to map the parameter space of the 
competition model to the lower dimensional parameter space of the predator- 
prey system. For simplicity we focus on a symmetrical situation in which both 
populations of predators behave identically. We set 

P = P = 0.5 

gi,x = 92, x = 9x ( 19 ) 

9l,y = 92,y = g y 
m l,y = ^2,y = "V 



While this situation is quite special, let us emphasize that it does not mean 
that we are talking about a single population of predators that we have split 
arbitrarily into two groups. Such a degenerate situation would, among other 
pathologies, lead to a rank deficient Jacobian. In truth, the two populations 
can be quite different. What Eq. (fl9|) implies is only that they have the same 
impact on the prey population and that they respond in the same way to 
perturbations. 
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Using this assumption the bifurcation condition, Eq. (JT81) simplifies to 



= s x (g y - m y )(g y - m y ) + 0.5g x m y (g y - m y ) + 0.5g x m y {g y - m y ) 
= s x (g y - m y ) 2 + g x m y (g y - m y ). 



Since g y ^ m y this is equivalent to 

= s x (g y - m y ) + g x m y (21) 

and therefore 

9x = s x - (22) 

m y 



This is exactly the condition for bifurcations of saddle-node-type that we found 
in the predator-system. Also, under the assumptions stated in Eq. ffTUj) . the 
condition for Hopf bifurcations in the competition model maps exactly onto 
the corresponding condition, Eq. (fTOj) . in the predator-prey system. 




1.0L0 

Fig. 2. Bifurcation diagram of a system of two predators competing for a common 
prey (see sketch in the top right corner). The generalized parameters describing 
predator 2 have been chosen to be slightly different (a y 2 = 1.050^1, m,2,y = 0.95mi )y ) 
52,x = 0.95g2,y) Predator 2 causes 55 % of the biomass loss of the prey. As in 
Fig. [pj the surfaces mark Hopf bifurcations points (red) and bifurcations points of 
saddle-node-type (blue) . Although the diagram is not identical to Fig. CD it shows a 
strong similarity. (Further parameters: s x = 0.5, g y = 1.) 

Let us now see what happens if we relax Eq. ( fl9i) . For this purpose we change 
all parameters describing one of the competitors slightly. The resulting bifur- 
cation diagram is shown in Fig. [2j While the bifurcation surfaces shift a little, 
the structure of the bifurcation diagram is qualitatively conserved: No new 
bifurcation surfaces appear or vanish in response to the change of parameters. 
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It is conceivable that the two diagrams in the figure could still be exactly iden- 
tical if we had chosen the mapping of parameters on the axes of the diagram 
more cleverly. This last point will be discussed in more detail in Sec. El 

To summarize: We have used the approach of generalized modeling to iden- 
tify a situation in which the local bifurcation diagrams of a predator-prey 
system is exactly identical to the corresponding diagram of a system describ- 
ing exploitative competition. While this situation, is admittedly special in the 
context of the generalized model, it can be found in the whole family of con- 
ventional models that satisfy Eq. ffl9l) . The equivalence of bifurcation diagrams 
implies that the dynamics of the two models is, at least locally, identical and 
therefore the coarse-grained (predator-prey) model and the full (exploitative 
competition) model have similar local dynamics. 

o o o op oop 

1 = \ f = f\ = txi = u/ = 

o o o o o o o 
o o o o o o?p 

♦ \ / t / \ w 

0=0 = 0=00= = 

♦ ♦ f\ \f JL 
o 0000 o 

Fig. 3. Symbolical representation of some examples of the equivalence relation be- 
tween bifurcation diagrams of different food webs. The sketches indicate food web 
topologies. If species on the same trophic level interact with the same neighbors 
and are described by the same generalized parameters they can be merged without 
altering the local bifurcation diagram. In the numerical verification of the relations 
shown here passive prey switching has been assumed. 



3 Equivalence of more complex food webs 



In the previous section we have shown that under certain conditions different 
systems may be described by exactly the same local bifurcation diagrams. This 
insight is not limited to the minimal examples studied in the previous section 
but can be observed in many food webs. The Jacobian of a very general food 
web was derived in Q and in a slightly different parameterization in jij. This 
Jacobian describes the local dynamics in complex food webs with nonlinear 
interactions and takes passive prey switching into account. We have compared 
the bifurcation diagrams of 43 different food webs numerically. Some of the 
results are shown symbolically in Fig. [3J 
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Based on our results we conjecture the following: If two or more populations 
in the system interact with exactly the same neighbors and their generalized 
parameters are identical in the sense of Eq. (1191) then the populations can be 
modeled by a single dynamical equation with the corresponding parameters. 
Again we emphasize that identical generalized parameters do not imply an 
identity of organisms. 

In the following we call two food webs equivalent if they satisfy the assump- 
tions of the conjecture. In this case we say that the food web with the higher 
number of populations is reducible to the food web with the lower number of 
populations. 




Fig. 4. Top row: Equivalence of complex food webs. If parameters are chosen ap- 
propriately (see text) the bifurcation diagrams of certain complex food webs are 
exactly identical. The respective topologies are indicated in the top right corner. 
Both diagrams are identical to the corresponding diagram of a four-trophic food 
chain. Bottom row: Example of two food webs that cannot be mapped to food 
chains in the proposed way. In all four diagrams the red, yellow, and green surfaces 
correspond to Hopf bifurcations and the blue surface to bifurcations of saddle-node 
type (Parameters: s x = 0.5, g^ y = 1, gi x = g x , Oi = r 3 where j is the trophic level 
of i, m^y = m y . For details see 0]) 
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In particular the conjecture implies that there is a class of food webs that are 
reducible to food chains of the same number of trophic levels. An example of 
such food webs is shown in Fig. [4] (top panels). If one of the assumptions of the 
conjecture is violated we generally find that the food webs are not equivalent 
under a mapping of the form of Eq. ffl9|) (see bottom panels of Fig. H|). In 
particular a food chain can not be reduced to another food chain with fewer 
trophic levels. 

Additionally, we observed that the bifurcation diagram of a given food web 
contains at least as many Hopf bifurcation surfaces as the food chain with 
the same number of trophic levels, while food webs that are not reducible to 
food chains generally have additional bifurcation surfaces. Two examples of 
such food webs that are not reducible to food chains are also shown in Fig. [4] 
(bottom panels). 

The examples discussed above show that the food web topology has a strong 
impact on reducibility. We therefore have to ask what happens if we alter 
the topology of a reducible food web. At which point will the transition to 
irreducibility occur? In our general food web model the topology is captured by 
a number of parameters, that encode the strength of topological connections 
in analogy to the parameter f3 in the competition model from the previous 
section. We can tune these parameters continuously between reducible and 
irreducible topologies. 

Let us illustrate this point in a simple example. We consider a system in which 
two competing populations of predators feed on a single prey population. Both 
predator populations are consumed by a population of top predators. In our 
model the prey species is population 1, the competing predators are popula- 
tion 2 and 3 and the top-predator is population 4 (see sketch in Fig. We 
assume that both populations of predators are described by similar general- 
ized parameters in the sense of Eq. (TTOT) . except that the biomass flow from 
population 3 to population 4 may be less than the biomass flow from popu- 
lation 2 to population 4. We can imagine that population 3 has some defense 
against the predator that population 2 lacks. We measure the efficiency of this 
defense, by a parameter 6 3j4 , which indicates the strength of the trophic link. 
If 6 3 4 = 1 the defense is ineffective, and individuals of population 3 £1X6 cLS 
likely to be consumed by the top predator as individuals of population 2. We 
assume that in this case natural mortality of individuals of population 2 and 3 
can be neglected. If we decrease 63^ the efficiency of the defense increases. In 
this case only a certain portion of the biomass loss of population 3 is caused 
by the top predator while the remainder is due to natural mortality. If 63^ = 
then the biomass loss of population 3 because of predation vanishes, and the 
entire biomass loss of population 3 occurs because of natural mortality. 

A bifurcation diagram of the system described above is shown in Fig. In this 
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Fig. 5. Top row: Comparison of the bifurcation diagram of a tri-trophic food chain 
(top right) with a food web of variable topology (top left). In the variable food web 
the strength of one trophic link is controlled by the parameter 63^. This link is shown 
in grey in the symbolic representation of the food web topology. For parameters see 
Fig. S] except for m y , which is 1.1 in the left diagram. Bottom row: sketches of 
two-parameter slices through the three-parameter diagrams. The two parameter 
diagram for the tri-trophic food chain at m y = 1.1 (A) is exactly identical to the 
diagram of the variable system for 63,4 = 1 (B). If 63,4 is decreased the bifurcation 
diagram changes (C) but remains qualitatively similar to A. At still lower values 
(D) the bifurcation diagram differs qualitatively from the diagram of the food chain 
for m y = 1.1 but is still similar to the food chain diagram at higher m y . Finally, at 
low values of 63,4 the bifurcation diagram of the variable system is no longer similar 
to the bifurcation diagram of the tri-trophic food chain for any set of parameter 
values (E). 



diagram we have chosen 63 4 as one of the axes. For 6 3 4 = the food web can 
be reduced to the tri-trophic food chain. Therefore, the parameter plane with 
63,4 = is exactly identical to a two-parameter bifurcation diagram of this 
food chain at the corresponding exponent of closure m y = 1.1. As we decrease 
63,4 the exact identity is lost but the topology of the bifurcation surfaces 
at first stays qualitatively similar. At a certain value of 63^ the bifurcation 
surface of saddle-node-type folds back and disappears. From this point the two- 
parameter bifurcation diagram that we find in the food web for a given value 
of 63,4 is no longer equivalent to the bifurcation diagram of the tri-trophic food 
chain at the corresponding value of m y . However, note that the two-parameter 
diagram is still qualitatively similar to a two-parameter diagram of the food 
chain at a higher exponent of closure m y . If 6 3 4 is decreased further then a new 
Hopf bifurcation surface appears. The bifurcation diagram is now qualitatively 
different from the bifurcation diagram of the food chain. In particular, double 
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Hopf-bifurcations, i.e. intersections of Hopf bifurcation surfaces as shown in 
Fig. [5] (E), can not appear in our tri-trophic food chain model at all. 

Above we have conjectured that local bifurcation diagrams of a class of food 
webs are exactly reducible to the corresponding diagrams of food chains with 
the same number of trophic levels. Figure [5] shows that in the example a larger 
class of food web appears to exhibit qualitatively similar local dynamics. These 
insights, if proved, can be used to extend certain results on food chains to 
classes of food webs. 



One question that has been discussed in the ecological literature for quite 
some time is whether chaotic dynamics occur in natural food webs 

@ ^Hill- 
On the one hand even some of the simplest ecological models exhibit chaotic 
dynamics (e.g. 12l;ll3l;ll4j| ) . On the other hand chaos has rarely been observed 
in nature 15; 16; 17; 18|. In the past it has been argued that chaotic dynamics 
could be detrimental for the survival of the participating populations and 
could therefore disappear in the course of natural selection. However, selection 
acts primarily on the level of the individual. And, even on the level of the 
population there is some evidence that chaos, if considered in a spatial or 
meta-population context, promotes persistence 19; 2(3; 21 1. Perhaps the most 
widely accepted view is now that chaos may be present in the dynamics but 



is difficult to detect because of environmental noise 22; 231. However, another 



opinion is that the complexity of large ecological systems somehow prevents 
chaotic dynamics [H; 0; 25 . 



Since the approach of generalized modeling focuses on local dynamics, global 
dynamics such as chaos is not directly accessible. However, we can draw con- 
clusions on certain features of global dynamics by considering bifurcations 
of higher codimension. Unlike Hopf and saddle-node bifurcations, which are 
of codimension 1, bifurcations of higher codimension have more than one bi- 
furcation condition. In comparison to codimension- 1 bifurcations the higher 
codimension bifurcations are often neglected, since the tuning of more than 
one parameter is necessary to find them in experiments or numerical investi- 
gations. 



An example of a codimension-2 bifurcation is the double-Hopf bifurcation, 
formed at the intersection of two Hopf bifurcation surfaces, which we have 
already encountered above. Mathematical investigations of this bifurcation 
have shown that there is generically a region in parameter space in which 



chaotic dynamics occur close to the double-Hopf bifurcation point [26J]. In 
this parameter region chaotic dynamics can therefore, at least transiently, be 
observed. 



In a previous paper we have shown that double Hopf bifurcations exist in food 
chains with more than three trophic levels. This implies that chaotic dynamics 
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can be found generically in these food chains at least in some parameter space. 
Our conjecture implies that a certain class of food webs is exactly equivalent to 
food chains of the same length. This means that these food webs also exhibit 
double-Hopf bifurcations and therefore chaotic parameter regions. Further- 
more, double-Hopf bifurcations should also exist in the much larger class of 
food webs in which the bifurcation diagrams are qualitatively similar to the 
bifurcation diagrams of food chains. Finally, irreducible food webs in which 
the local dynamics differs qualitatively from that of food chains have, in our 
experience, a more complex bifurcation structure than food chains of compa- 
rable length and therefore the corresponding bifurcation diagrams also contain 
double-Hopf bifurcations. 

From the evidence given above we conclude that chaotic parameter regions 
should exist in a large class of complex food webs with more than three trophic 
levels. It is therefore unlikely that the complexity of natural food webs in itself 
leads to the avoidance of chaotic dynamics. 



4 Remarks and Discussion 

In this paper we have used the approach of generalized models to show that 
in certain food webs the local bifurcation diagrams are exactly identical to 
the corresponding bifurcation diagrams in other food webs and food chains. 
Moreover, we have conjectured that there is a universal rule that governs 
which food webs show this equivalence. This conjecture, if proved, can be used 
to generalize certain findings on food chains to classes of food webs. Perhaps 
more importantly, it holds the promise to yield an analytical procedure for the 
reduction of complex food-web models or even dynamical biological networks 
in general. 

A method for coarse-graining food webs based on the equivalence of local 
bifurcations had several advantages: It would probably preserve local dynamics 
exactly, and also preserve certain features of the global dynamics through 
conservation of local bifurcations of higher codimension. Moreover, it would 
probably provide hard criteria governing which populations could be coarse- 
grained into a single variable and which could not. 

Before the present observations can be extended to a method for model reduc- 
tion more research is certainly necessary. This work should aim to understand 
the equivalence, prove the stated conjecture and in particular explore its va- 
lidity in the large class of systems in which we so far only observe a qualitative 
equivalence of bifurcation diagrams. Let us therefore conclude with some re- 
marks that we believe will be conductive for this effort. 
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One feature of the mapping of a food web to a simpler food web is that it 
reduces the dimensionality of both the state space and the parameter space. 
In the context of our approach the mapping in state space is not of much 
concern as we do not need to specify it explicitly. By contrast the mapping in 
parameter space is more intriguing; We have to specify how the parameters 
of the original model are translated to the lower number of parameters in the 
reduced model. 

In this paper we have postulated the mapping of parameters ad hoc based on 
biological reasoning and only subsequently verified that it produced the desired 
result. For instance in our discussion of the system of competing predators we 
have demanded that both predators exhibit the same exponent of closure 
m y = m l y = m 2) y (Eq. [T^]) . Only thereafter we could compare the bifurcation 
diagram to that of the predator-prey system. For future investigations a dif- 
ferent approach seems more promising: Instead of demanding the identity of 
certain parameters and then checking for the identity of bifurcation diagrams, 
we could ask which mapping of parameters, say m y (mi iy , m 2 , y ) will produce 
a bifurcation diagram that is identical to the one of the simpler system. Of 
course we do not need to restrict this mapping to parameters between which 
a direct biological correspondence exists, but could allow the parameter of 
the reduced system m y also to depend on other parameters such as gi tY and 
<?2,y While more mathematical work is certainly necessary, we believe that it 
should be possible to derive a rigorous procedure for the parameter mapping. 

We believe that the procedure described above will reveal exact reductions 
for the large class of systems in which we so far only observe a qualitative 
similarity of bifurcation diagrams. Nevertheless, it is obvious that for certain 
pairs of food webs there cannot be any mapping of parameters that leads to 
identical bifurcation diagrams. For instance, we have seen that the bifurcation 
diagram of the four-trophic food chain contains a double- Hopf bifurcation line. 
Since this bifurcation is characterized by four purely imaginary eigenvalues of 
the Jacobian it simply can not be accommodated in any system with less than 
four state variables. Therefore the local bifurcation diagram of the four-trophic 
food chain can never be mapped on that of the tri-trophic food chain. 

The fact that generalized modeling reveals bifurcations of higher codimension 
with relative ease is advantageous as it provides examples of these bifurcations 
for mathematical analysis and allows us to gain insights into global features 
of the dynamics based on a local analysis. In the present paper this has be 
utilized to provide evidence for chaos in a large class of food webs. This and 
other bifurcations also reveal the presense of global homoclinic and heteroclinic 
bifurcations. It is conceivable that future progress in normal form theory will 
enable us to extract even more information from the bifurcations of higher 
codimension. 
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So far we have focused on generalized models while in the literature most mod- 
eling is done by conventional models in which the interactions are restricted to 
specific functional forms. Note however, that it is straight forward to step back 
and forth between generalized and conventional representations of a system. 
By replacing the specific functional forms by general functions, a conventional 
model can be turned into a generalized model and vice versa 23]. The corre- 
sponding mapping between generalized and specific model parameters can be 
computed analytically and thus provides a direct relation between the model- 
ing approaches [28]. For a smoother transition, it is possible to study hybrid 



models in which some functional relationships are specified while others are 
kept general 



Let us emphasize that the identity of bifurcation diagrams is therefore not 
restricted to the realm of generalized models; One should be able to observe 
the same identity also in conventional models. However, in conventional mod- 
els additional difficulties arise since the mapping is affected by the position 
of steady states in state space which is often prohibitively difficult to study 
analytically. In the past we have pointed out that the parameters that are 
identified in generalized modeling are in a certain sense more natural than the 
parameters that are commonly introduced in conventional modeling 0;0]. In 
the present paper this enabled us to spot the identity of bifurcation diagrams 
even based on very simple ad hoc assumptions. While the same identity should 
exist also in conventional models the corresponding mapping of the conven- 
tional parameters is probably much more involved and hence the identity is 
more difficult to spot. 



To conclude: We have shown that in generalized models an exact equivalence 
between the bifurcation diagrams of different food webs can be observed, albeit 
in a relatively special case. Several results seem to indicate that this identity 
can be extended to a much larger class of systems if more refined parameter 
mappings are used. We hope that this insight will eventually evolve into a 
method for model reduction. In the context of such a method an intermediate 
coarse-graining step with generalized models will probably be useful even if 
the desired result and starting point are conventional models. 
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